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^ . interface tension on the coarse graining scale k. A stable range appears only if the 

renormalized dimensionless couplings at the critical temperature are small. This 
gives a quantitative criterion for the validity of computations within Langer's theory 
of spontaneous bubble nucleation. 
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The discussion of the dynamics of a first order phase transition [lj usually relies on 
the study of a non-convex potential or free energy. The decay of unstable minima is 
associated either with tunneling fluctuations through barriers in the potential 0, or, at 
non-zero temperature, with thermal fluctuations above them ||. However, the effective 
potential (4|, which seems at first sight a natural tool for such studies, is expected to be 
a convex quantity with no barrier. The resolution of this paradox lies in the realization 
that the effective potential is convex because the tunneling or thermal fluctuations are 
incorporated in it. These fluctuations are associated with low frequency modes, while the 
non-convex part of the potential is related to the classical potential and the integration 
of high frequency modes. A natural approach to the study of first order phase transitions 
separates the problem in two parts. First, the high frequency modes are integrated out, 
with the possible generation of new minima through radiative symmetry breaking [[J. 
Subsequently, the decay of unstable minima is discussed with semiclassical techniques 
[0, [3], in the non-convex potential that has resulted from the first step. This leads us to 
the notion of the coarse grained free energy, which is fundamental in statistical physics. 
Every physical system has a characteristic length scale associated with it. The dynamics 
of smaller length scales is integrated out, and is incorporated in the parameters of the 
free energy one uses for the study of the behavior at larger length scales. 

The notion of coarse graining is absent in the perturbative approach to the calculation 
of the effective potential ||. This is the main reason for the non- convergence of the 
perturbative series near the maxima of the classical potential, and the appearance of 
imaginary parts in the perturbative effective potential. Despite attempts to give a physical 
interpretation to these imaginary parts ||, a satisfactory discussion of tunneling must 
incorporate the notion of coarse graining. The Wilson approach to the renormalization 
group provides the appropriate framework |7| . We employ here the method of the effective 
average action I\. , which results from the integration of fluctuations with characteristic 
momenta larger than a given scale k. The dependence of on k is described by an exact 
renormalization group equation^ ||, ^] . For large values of k (of the order of the ultraviolet 
cutoff A of the theory) the effective average action is equal to the classical action (no 
fluctuations are integrated out), while for k — > it becomes the standard effective action 
(all fluctuations are integrated out). For non- zero k the effective average action has the 
properties of a coarse grained free energy. Its non-derivative part (the effective average 
potential Z7&) is not necessarily convex. The coarse graining scale can be identified with 
k. 

In this letter we provide an explicit demonstration of how such a potential can be 
obtained starting from the microscopic or classical action of a field theory. We investigate 
the dependence of the effective average potential t/& and the 'classical' surface tension a k 
on the coarse graining scale k with special emphasis on the question of the validity of 
Langer's theory of bubble formation. We study the first order phase transitions for the 
Abelian Higgs model and for a scalar matrix model in three dimensions. An application 
of the formalism to the case of the high temperature phase transitions for the Abelian and 

1 For related work see |@, ffij fy) . 
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SU(2) Higgs models is given in ref. [12[| and a discussion of the first order phase transition 



in matrix models can be found in ref. 13 



Near a phase transition, the three dimensional Abelian Higgs model describes the 



behavior of ordinary superconductors flTJj]. It can also be viewed as the effective theory 
resulting from the non-zero temperature four dimensional model near the critical temper- 
ature. The dependence of the effective average potential U k {p) and the running renor- 
malized gauge coupling e^k) on the coarse graining scale k is governed by the evolution 
equations M 



dU k {p) 
dt 



de 2 R ( k ) 
dt 



1 



d 3 g 8P k 

2(2tt)3 dt {P k (q) + U' k (p) + 2U'M P 

1 2 \ 

+ 



Pk(q) + U' k (p) P k (q) + 2e 2 R (k) P/ 
0.84 e%(k) 



6tt 5 



where t = \n(k/A) (with A the ultraviolet cutoff of the theory) and p = \<fi\ 2 /2 (with the 
complex order parameter). Primes denote derivatives with respect to p. The three terms 
on the r.h.s. correspond to the contributions of the radial and Goldstone scalar modes 
and the gauge field. The inverse propagator 



Pk{q) 



exp (—q 2 /k 2 ) 



provides for an infrared cutoff which acts like a mass term 
q 2 <C k 2 , while it leaves unaffected the modes with q 2 ^> k 2 . 



(2) 

^ k 2 for the modes with 
The momentum integral 



on the r.h.s. of eq. (|l|) can be written in terms of dimensionless functions Iq{w), whose 
arguments are given by the rescaled mass terms (U' k (p) + 2U k (p)p)/k 2 , U' k (p)/k 2 and 
2e 2 R (k)p/ k 2 . These functions fall off for large values of w, following a power law. As a result 
they introduce threshold behavior, which leads to the decoupling of massive modes from 
the evolution equations || |TJ, 16]. The derivation of eq. ([I]) under some approximations, 
starting from the exact renormalization group equation for the effective average action, 
is given in ref. |], [TJ], 0. The approximations concern the omission of the anomalous 
dimension of the scalar field, the effective field dependence of the gauge coupling and the 
higher derivative terms in the action. The evolution starts at k = A, where the effective 
average potential is equal to the microscopic or classical one U^(p) = o^a(p ~~ Poa) 2 
and the running gauge coupling is equal to the bare coupling e\. In the opposite limit 
k — ► 0, Uk(p) becomes equal to the (convex) effective potential U(p) = Uo(p) and the 
gauge coupling assumes its renormalized value en = eu(k = 0). Two algorithms for the 
numerical integration of eq. (|]) have been presented in detail ||17|| . The phase transition 
is approached by fixing Aa and e\ and tuning p A- The system exhibits a second order 
phase transition for = which corresponds to the Wilson-Fisher fixed point of the 
0(2) symmetric Heisenberg model. For large enough e\/\^ the phase transition is first 
order M, HI O 
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In fig. 1 we display the solution of eq. (U) for Aa = 0.01A, e\ = 0.1A and poA — 0.867A. 
All the quantities in the figures are expressed in units of the ultraviolet cutoff A. Initially 
the potential has only one minimum away from the origin, which evolves proportionally 
to the coarse graining scale k. At some point a new minimum appears at the origin. It 
is induced by the integration of fluctuations, through the generalization of the Coleman- 
Weinberg mechanism. The evolution slows down at the later stages, and for k/A around 
0.02 the potential converges towards a stable non-convex profile with two minima of 
equal depth. Around the minima the scale k becomes smaller than the mass of the various 
massive modes, and this induces their decoupling. We have stopped the evolution at a non- 
zero kf, for which the shape of the potential near the minima is stable. The presence of the 
non-convex part is explained by this non-zero value of k. We have not yet integrated out 
all the fluctuations, which should render the effective potential convex. More specifically, 
the fluctuations which interpolate between the two minima of fig. 1 are not included 
effectively in the non-convex potential. They are the ones that trigger the tunneling and 
drive the first order phase transition. If we continue the evolution all the way to k = 0, 
these interpolating configurations will be gradually integrated out. As a result, the height 
of the barrier will start getting smaller, until the region of the potential between the two 
minima becomes flat \TE\. The evolution of the characteristics of the potential is depicted 
in fig. 2. We plot the location of the minimum away from the origin p m in, the value of the 
potential at the minimum {Uk)min, the location of the maximum p max , the value of the 
potential (U k ) max and the curvature (d 2 U k /d(p 2 ) max = U' k (p max ) + 2p max U' k r (p max ) at the 
maximum. We observe that these parameters have almost constant values in the region 
k/A ~ 0.02 - 0.03. 

In fig. 3 and 4 we present the effective average potential and its characteristics for 
parameters corresponding to a weaker first order transition. For Aa = 0.1 A, e\ = 0.1 A 
and Poa — 0.1 71 A the discontinuity in the scalar field expectation value is about nine times 
smaller than for fig. 1. The most important difference is that the potential never becomes 
relatively stable for a range of k. During the later stages of the evolution, its outer part 
(for scalar field values larger than the location of the minimum) starts approaching a 
stable profile, due to the decoupling of the massive modes in this region. However, in the 
same range of k the non-convex part starts already becoming flatter, as configurations 
interpolating between the two minima are being integrated out. The negative curvature 
at the top of the barrier is expected |15| to behave ~ —k 2 during this stage. This has 
been verified explicitly through the analytical integration of the evolution equation for the 



O(N) symmetric scalar theory in the large N limit [19 . We clearly observe the onset of 
this behavior of d 2 U /dip 2 in fig. 4 in a range of k before p m i n settles. Also the maximum of 
Uk decreases before the minimum settles. In this case it is far from obvious which coarse 
graining scale k should be chosen for a definition of important nucleation characteristics 
such as the interface tension. 

The behavior of the coarse grained effective potential for a fluctuation induced first 
order phase transition, as considered above, is not particular to the Abelian Higgs model. 
For example, our discussion can be extended with minor modifications to the electroweak 
phase transition for the range of Higgs-scalar masses where it is first order flIS|. It also 
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can be observed in pure scalar theories. One may consider models with U(N) x U(N) 
symmetry with a scalar field in the (N, N) representation, described by a complex N x N 



matrix [13]. The cases N = 2, 3 have an interesting relation to high temperature 



strong interaction physics and chiral symmetry breaking |20| and non-perturbative flow 
equations have been studied in this context We will concentrate here on N = 2. The 
most general effective average potential Uk(p,r) can then be expressed as a function of 
only two invariants, namely 

p = tr(0V), r = 2tr(Vv-^p) • (3) 

The microscopic or classical potential U\ for these models can be characterized by two 
quartic couplings Aia, A 2 a and a mass term {jx\ > 0), 

X _ i_ 

U A (p, t) = -fl 2 A p + - Aiap 2 + -A 2A r . (4) 

In the limit Aia ~~ ► oo, A 2 a ~~ ► °o this also covers the model of unitary matrices. One 
observes two symmetry breaking patterns for A 2 a > and A 2 a < respectively. The case 
A 2 a = denotes the boundary between the two phases. In this special case the theory 
exhibits an enhanced 0(8) symmetry and one finds a second order phase transition. For 
the symmetry breaking pattern £7(2) x £7(2) — > U{2) (A 2 a > 0) the phase transition 
is always first order. In this case the relevant information for the phase transition is 
contained in t4(p) = Uk(p,r = 0). The discussion of the dependence of the effective 
average potential Uk(p) on the coarse graining scale k can be presented along the same 
lines as for the Abelian Higgs model and the relevant flow equations can be found in 
ref. |TJ|]. Here the second quartic coupling A 2 a for the scalar model plays the role of the 
gauge coupling e\ in the Abelian Higgs model. In addition to the above treatment of the 
Abelian Higgs model, the employed approximation for the scalar model takes into account 
a ^-dependent wave function renormalization constant Zk for the fields and the effective 
field dependence of the second quartic coupling. 

In the following we will use the scalar matrix model to establish a quantitative criterion 
for the validity of the standard treatment of spontaneous bubble nucleation as described 
by Langer's theory ||T|. Langer's approach relies, on the one hand, on the definition 
of a suitable coarse grained free energy 1^ with a coarse graining scale k and, on the 
other hand, on a saddle point approximation for the treatment of fluctuations around 
the 'critical bubble'. The problem is therefore separated in two parts: The first part 
concerns the treatment of fluctuations with momenta q 2 ^ k 2 which are included in the 
coarse grained free energy. The second part deals with an estimate of fluctuations around 
the bubble, for which only fluctuations with momenta smaller than k must be considered. 
To be explicit we consider a spherical bubble where the bubble wall with thickness A is 
thin as compared to the bubble radius R, i.e. A <C R. In this thin wall approximation 
the bubble nucleation rate f, which describes the probability per unit volume per unit 
time for the transition to the new vacuum, can be written in the form || [| 

r = Aexp(-M. (5) 
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Here the 'classical' surface tension a k is given in our conventions by 



a k = 2 d^J2Z k U k {ip) (6) 



o 



where ip = (p/2) 1 / 2 and (p denotes the second zero of U k (<p) near the outer minimum at 
fmin- For the difference in the free energy density e we include fluctuations with arbitrarily 
small momenta, 

e = lim(£/fc(0) - U k {ip min )). (7) 
fc— >o 

In contrast, the long wavelength contributions to the true surface tensionP] are effectively 
cut off by the characteristic length scale of the bubble surface. For the free energy of the 
critical bubble the modes with q 2 <; k 2 are incorporated in 1^ = 167Rx|/3e 2 (lowest order 
or classical contribution) and they influence a k . The modes with q 2 ^ k 2 contribute to the 
'fluctuation determinant' A k , which also contains dynamical factors. Here A k depends 
on k through the effective ultraviolet cutoff for these fluctuations, which is present since 
fluctuations with momenta larger than k are already included in T k °\ A more general 
discussion which does not rely on the thin wall approximation or a saddle point approx- 
imation is given in ref. fll3"| . Langer's formula for bubble nucleation amounts essentially 



to a perturbative one loop estimate of A k . It is clear that k is only a technical construct 
and for physical quantities like the bubble nucleation rate the fc-dependence of In A k and 
Ik must cancel. A strong dependence of a k on the coarse graining scale k is only com- 
patible with a large contribution from the higher orders of the saddle point expansion. 
The fc-dependence of o k therefore gives direct information about the convergence of the 
saddle point approximation and the validity of Langer's formula. 

We consider in detail the dependence of the surface tension a k on the coarse graining 
scale k at the phase transition (e = 0) for three examples. They are distinguished by 
different choices for the quartic couplings Aia and A 2 a of the short distance potential U\ 
given by eq. (|J). The choice Aia/A = 0.1, A 2 a/A = 2 corresponds to a strong first order 
phase transition with renormalized masses not much smaller than the cutoff scale A. In 
this case the fc-dependence of the effective average potential resembles the one presented 
in fig. 1. In contrast we give two examples where the dependence of a k on the coarse 
graining scale becomes of crucial importance. The choice Aia/ A = 2, A 2 a/A = 0.1 leads 
to a weak first order phase transition with small renormalized masses and the behavior 
of the effective average potential is similar to the one given in fig. 3. The coarse grained 
potential and the surface tension show a high sensitivity on the scale k. A more increased 
sensitivity on the scale k can be observed for Aia/A = 4, A 2 a/A = 70 which corresponds 
to a relatively strong first order phase transition. 

The ^-dependence of the surface tension a k is displayed in fig. 5. Here a k is normal- 
ized to its maximum value cr* where <r*/A 2 = 1.67 x 10 _2 (8.41 x 10~ n )(1.01 x 10~ 3 ) 
for Aia/A = 0.1(2)(4), A 2A /A = 2(0.1)(70). The scale k f is given by \m max \ with 

2 The true surface tension is a 'measurable quantity'. It is independent of k and all fluctuations must 
be included. It therefore differs, in general, from a k which includes only part of the fluctuations. 
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m max = Z^k 1 Pmax(d 2 Uk/dp 2 )(p max ) denoting the renormalized mass term at the top 
of the potential barrier p max . More precisely, we chose (kj — \'m 2 nax {kf)\)/k 2 f = 0.01. For 
Aia/A = 0.1, A2A/A = 2 the curve exhibits a small /c-dependence around its maximum 
and cr* ~ (Tfe.. For the second and the third example one observes that o kf becomes 
considerably smaller than the maximum value due to a strong fc-dependence. The coarse 
graining scale k should not be taken smaller than the inverse bubble wall thickness A -1 
[jnj. This ensures that the detailed properties of the bubble are irrelevant for the com- 
putation of the 'classical' surface tension 07.. We estimate the bubble wall thickness A 
by 

A = 2Z k ^ (8) 
o~k 

where we have taken the gradient energy as half the total surface energy and approximated 
the mean field gradient at the bubble wall by (p min /A. For the given examples we observe 
A^(kf) ~ kf/A. We choose k ~ kf as the coarse graining scale. 

In order to quantify the differences between the three examples we have displayed 
some characteristic quantities in the table. The renormalized couplings 



A lij = Z kf {Pmin) , A 2 ft = 4Z fc/ — {Pmin) (9 J 

are normalized with respect to the mass term 

m R = {2Z kPmin \ 1R ) 1 ' 2 . (10) 
In addition we give the mass term 

m 2R = {ZkPmin^R) 1 ^ 2 (H) 

corresponding to the curvature of the potential in the direction of the second invariant 
t. All couplings and masses are evaluated at the critical temperature (critical P o\). In 
comparison with fig. 5 one observes in the vicinity of kf a weaker scale dependence of a k 
for smaller effective couplings. In particular, a reasonably weak scale dependence of 
and <Tfc requires 

^ « 1 • (12) 
m R 

This establishes a quantitative criterion for the range where Langer's theory can be used 
without paying too much attention to the precise definition of the coarse graining. In the 
table we also present kf/A and the renormalized masses in units of A, which indicate the 
strength of the phase transition. In particular, for the relatively strong phase transition of 
the third example with slightly larger effective couplings one observes an increased scale 
dependence as compared to the weak phase transition of the second example. This clearly 
shows that the 'strength' of the phase transition is, in general, not the primary criterion 
for the applicability of Langer's theory. 

In addition to the dependence on k, the coarse grained free energy depends also on 
the precise shape of the infrared cutoff function or the inverse average propagator Pk{q) 
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given by eq. (0). Analytical studies in the Abelian Higgs model indicate [22 that this 
scheme dependence is rather weak for the effective potential and the surface tension. 

In summary, we have shown that the coarse grained free energy cannot be defined 
without detailed information on the coarse graining scale k unless the effective dimen- 
sionless couplings are small at the phase transition. For Abelian and non- Abelian gauge 
theories with a small gauge coupling at the scale k = A this coincides with a relatively 
strong first order transition. Only for small couplings a range with a weak fc-dependence 
of the classical surface tension appears. There is a close relation between the depen- 
dence of the coarse grained free energy on the coarse graining scale and the reliability 
of the saddle point approximation in Langer's theory of bubble nucleation. For a strong 
fc-dependence of U k a small variation in the coarse graining scale can induce large changes 
in the predicted nucleation rate in lowest order in a saddle point approximation. In this 
case the fc-dependence of the prefactor A k has also to be computed. Therefore, for strong 
dimensionless couplings a realistic estimate of the nucleation rate needs the capability 
to compute InA^ with the same accuracy as 16iral/3e 2 and a check of the cancelation 
of the /^-dependence in the combined expression Our observation that the details of 
the coarse graining prescription become less important in the case of small dimensionless 
couplings is consistent with the fact that typically small couplings are needed for a reliable 
saddle point approximation for A k . 
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Figures 

1. Evolution of the potential for a strongly first order phase transition. 



2. Characteristics of the potential for a strongly first order phase transition in depen- 
dence on the coarse graining scale k. 

3. Evolution of the potential for a weakly first order phase transition. 

4. Characteristics of the potential for a weakly first order phase transition in depen- 
dence on the coarse graining scale k. 

5. The normalized surface tension a k /a^ as a function of \n{k/kf). The short distance 
parameters are (1) A iA /A = 0.1, A 2A /A = 2, (2) A iA /A = 2, A 2A /A = 0.1, (3) 
Aia/A = 4, A 2A /A = 70. 



Table 

1. Effective dimensionless renormalized couplings \\R/m R and A 2 j?/m fi at the phase 
transition. The critical couplings and mass terms m c R , m\ R are evaluated at the 
scale kf. 
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